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Abstract 

We consider the energy weighted sums (EWS) for isovector giant dipole resonances (IVGDR) 
in finite nuclei within Landau kinetic theory. The dependence of both IVGDR energy, -EtvgdRj 
and the EWS enhancement factor, k(A), on the mass number A occurs because of the boundary 
condition on the moving nuclear surface. The values of -Ejvgdr^ 1 ^ 3 and n(A) increase with A. The 
obtained value of the enhancement factor is about 10% for light nuclei and reaches approximately 
20% for heavy nuclei. A fit of the enhancement factor to the proper experimental data provides a 
value for the isovector Landau amplitude of F[ ~ 1.1. 

PACS numbers: 21.60.Ev, 24.30.Cz 
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I. INTRODUCTION 



A microscopic description of the isovector excitations within the nuclear Fermi-liquid 
theory requires the use of two Landau amplitudes Fq and F[ to model the nucleon-nucleon 
interaction. The amplitude Fq provides the isospin symmetry energy, whereas the inclusion 
of the velocity-dependent force ~ F[ leads to the renormalization of the isovector energy 
weighted sums (EWS) jl, 2]. In particular, the isovector EWS m\ becomes dependent 
on the nucleon-nucleon interaction. This is opposite to the case of the isoscalar excitations 
where the corresponding sum wi\ is model independent. Moreover, the presence of the 
velocity dependent force gives rise to an exceed of the 100% exhaustion of Thomas-Reiche- 
Khun (TRK) sum rule for the isovector giant dipole resonances (IVGDR). The origin of the 
corresponding enhancement factor of the IVGDR sum mi was intensively investigated for 
both the nuclear matter and the finite nuclei, see Refs. Q, Q] and references therein. The 
RPA calculations of the enhancement factor for the symmetric nuclear matter were recently 
performed for a representative set of Skyrme forces in Ref. (J). As shown in 4], the value of 
the enhancement factor is changed in almost 2 times depending on the choice of the Skyrme 
force parametrization. The high sensitivity of_the enhancement factor to the choice of the 
Skyrme forces was also demonstrated in Ref. 



5| for the nucleus 208 Pb. 



In finite nuclei, both the IVGDR eigenenergy -Eivgdr and the EWS m,\ are rather com- 
plicated functions of the mass number A. In contrast to the classical Steinwedel- Jensen 
model 6], the value -EWgdr • A 1 ^ 3 is not a constant but increases with A. The theoretical 
approaches to microscopic description of IVGDR are mainly based on the Random Phase 
Approximation (RPA) [7|. The RPA analysis of the enhancement factor for some spherical 
nuclei (but not its A- dependence) has been recently done in Ref. js] 

overestimate the magnitude of the enhancement factor, 



Note, however, that 

the RPA calculations 

see below FIG. [3j Note also, that within the RPA, the highly excited IVGDR is strongly 
fragmented over a wide energy interval and a special averaging procedure has to be applied 
to pick up a smooth A-dependence of the IVGDR characteristics. Furthermore, the RPA 
calculations of the EWS are restricted due to taking into account lp — lh excitations 
only and, thereby, one can expect an underestimation of the contribution to mi from more 



complicated states 



id, 



In this work, we study the IVGDR within the Landau kinetic theory [121 . that is extended 
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to the finite two-component Fermi-liquid drop. The A-dependence of both the IVGDR 
eigenenergy and the corresponding EWS occurs due to the boundary conditions on the 



moving nuclear surface. Our approach is more general than the scaling model [13j or the fluid 



dynamic approaches 14[ due to the fact that it takes into consideration all multipolarities 
of the Fermi surface distortions. 

In Sect. II, we start from the response theory based on the collisionless kinetic Landau- 
Vlasov approach to the isovector excitations. We derive the main characteristics of the 
IVGDR, taking into account the velocity dependent part of the Landau isovector interaction. 
In Sect. Ill, we obtain the boundary condition for the isovector sound mode on the free 
nuclear surface. In Sect. IV, the numerical calculations for the IVGDR eigenenergy and the 
corresponding EWS sum m,\ are presented. Conclusions are given in Sect. V. 

II. RESPONSE FUNCTION AND ENERGY- WEIGHTED SUMS 

We consider the response of the nucleus to an external periodic in time field U ext (t): 

U cxt (t)=X e- lujt q + X*e^q*, (1) 

where q is the Hermitian one-particle operator, which depends on both spatial and isospin 
coordinates. If we assume Aq <C 1, the quantum mechanical expectation of the operator q is 



given by the following form 15] 



(q) = X (uj)Xoe~^ + X*(u J )X*e^, (2) 

where x(^) is the linear response function. We will evaluate the isovector density-density 
response function x{ u ) by using 

A 
3=1 

where Tj = 1 for the neutron and Tj — — 1 for the proton. 

To evaluate the response function x( u )i we use the collisionless kinetic Landau- Vlasov 
equation for the isovector excitations 

^Sf + v- V r 5f - V p / cq • V r (6U scU + U cxt ) = 0. (3) 

Here, 5f is the isovector variation of the Wigner distribution function, v is the nucleon veloc- 
ity, /eq is the equilibrium distribution function and SU se \{ describes the dynamical component 
of the selfconsistent mean field. 
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The solution to Eq. ([3]) can be written in terms of a plane wave [lz 

5f = -^u g ip)e^ t \ (4) 



where e p = p 2 /2m*, m* is the nucleon effective mass and u^(p) is an unknown amplitude. 
The variation of the self-consistent field 5U se \{ in Eq. (jHJ) is obtained through the isovector 
interaction amplitude F'(p,j?) as 



5U< 



self 



(5) 



The interaction amplitude F'(p,p^) is usually parametrized in terms of the Landau constants 

Q 

j oo 

F\p^) = —^F' k P k (p-p), p = p/p, 



F' k as 



(6) 



k=0 



where Pk{x) is the Legendre polynomial, Np is the averaged density of states at the Fermi 
surface, given by 

N F = ~[^L^=^ (7) 



J (27rfr) 3 de p tt 2 H 3 
and Pf is the Fermi momentum. 

The parameter Fq is related to the isotopic symmetry energy C sym in the Weizsacker mass 
formula jf], Q 



C sym =-e F (l + F^), (8) 
with the Fermi energy Bp — p 2 F /2m*. In the following, we will assume that 



f^o, f{^o, f; >2 = o. 



(9) 



By using Eq. (j3J) to get the amplitudes v$(p), we obtain the isovector response function 
in the form 17 1 

Qoo(s) 



i-g(s) Qoo( s ) 



Here, s = um*/qpp, 



Qoo( s ) = n fQoo(s), Qoo(s) = 1 + - In 



s - 1 



s+1 



.7T 



+<- a e(i-H), 



and 



(10) 



(11) 



(12) 
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The frequencies of isovector eigenvibrations [the poles of the response function (flQl) ] can 
be derived from the dispersion relation 



1 - g(s) Q 00 (s) = 0. (13) 

The response function in Eq. fflOl) allows us to evaluate the isovector EWS 

1 f°° 

m k = - d(Huj) (hu) k lmx(uj). (14) 
7T Jo 

Using the dispersion relation between lmx{uj) and Rex(w), and the asymptotic behavior of 
Rex(^) at <jj — > and — ► oo limits, one obtains (see also Ref. [l|) 

m - 1 = ^^ mi = ^ 2 ' m3 = ^^- (15) 

Here, we introduced the renormalized (because of the Fermi surface distortion effect) isotopic 
symmetry energy C' sym = C sym + 85^/15 and the effective mass m! = m/(l + knm) for the 
isovector channel, where k^m is the enhancement factor of the sum rule (for nuclear matter), 
which is defined by the relation 

1 + k nm = (1 + F{/3)/(l + F 1 /3), (16) 

where is the Landau amplitude for the isoscalar channel. In contrast to the isoscalar 
excitations, the isovector EWS sum mi in Eq. ffloT) is not model independent in the sense 
that it depends on the effective mass m' and, thereby, on the interaction amplitudes F\ and 
F[. We also point out that the EWS of Eq. ffTBT) allows us to evaluate the constrained 
energy, E constT , and the scaling energy, E sc , for the IVGDR. Namely, 



w m_! V ml y mi V m 



III. BOUNDARY CONDITION 

For finite nuclei, the dispersion relation [Eq. fll3p ] has to be supplemented by a corre- 
sponding boundary condition. The boundary condition can be viewed as a condition for a 
balance of the forces acting on the free nuclear surface 



n ■ F 



+ n-F s = 0, (18) 

s 
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where n is the unit vector normal to the nuclear surface S, the internal force F is associated 
with the isovector sound wave and F$ is the isovector surface tension force. The internal 
force F is defined by the momentum fli 
F a = np5P a/ 3 where bP a p is given by 



ux tensor SP a 3 inside the nuclear volume. Thus, 



SR 



a/3 



P-F (V a £/3 + V^ Q ) + ( C sym p eq - -p F ) V • £ 6 a p. 



(19) 



Here, £ is the displacement field, p eq = (p n ,cq + Pp,eq)/2, and 



» F ~ 2 Peq£F TTi^3 



3s 2 



(20) 



The surface force Fs in Eq. ( |T8|) is defined as F, j; 



n,j 8P^ where is the pressure 
caused by the isovector polarizations at the nuclear surface [13J, ll7| , given by 



SPs 



O Peg 

3 r 



Q 5R X 



(21) 



Here, ro is the mean distance between nucleons, Q is related to the surface symmetry energy 
in the Weizsacker mass formula and 5R\ = Roats(t)Yi (r). The amplitude as(t) of the 
isovector vibrations of the nuclear surface is connected to the corresponding amplitude £ of 
the displacement field in a sound wave. 

Finally, from Eqs. (JT51) . (1191) and (1211) we derive the following secular equation for the 
wave number q 

2p eq Q 



Peq r PF PF 

4 sym+ 3 x 2 . 



Pf _ 
x 3gr (l + «atm) 



0. 



(22) 



where x = qRo and i?o — r^A 1 ^. In the limit Q — > oo, the boundary condition [Eq. (|22l) ] 
gives rise to the boundary condition j{(x) = of the Steinwedel- Jensen model 6|. We point 
out that the secular equation (122!) for g has to be solved consistently with the dispersion 
equation ( 1T3|) for s. 



IV. NUMERICAL CALCULATIONS 



Numerical calculations were carried out by using the following set of nuclear parameters: 
ro = 1.2 fm, Fi = —0.3 and Fq = 1.41. According to Eq. ([8]), the bulk symmetry energy 
C sjm is equal to 60 MeV [6j. The value of the Landau amplitude F[ will be discussed in the 
following. 
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In FIG. [U the solid curve 1 shows the A dependence of the value x obtained from the 
secular equation (1221) ; it is consistent with the dispersion relation [Eq. ([TBI ] including all 
multipolarities / of the Fermi surface distortions. Curve 3 represents the analogous result 
but for the velocity- independent nuclear forces (i.e., for Fx — and F[ = 0). Curve 2 
shows a solution to the secular equation (1221) . when one takes into account the Fermi surface 

n 

distortions up to quadrupole order (scaling approximation with I < 2) [13j. In the last case, 
instead of solving the dispersion equation ffTB"]) . we have used the expression s 2 = (9/5+Fq)/3, 
which gives the dimensionless sound velocity s for I < 2 As seen from FIG. [I], the value 
of x is significantly smaller than the corresponding one obtained within the Steinwedel- 
Jensen model (dashed line in FIG. [T]). 

In FIG. 12] we show a dependence of the IVGDR energy (multiplied by the factor A 1 / 3 ) on 
the mass number A. The calculations have been performed at Q = 10.5 MeV and F[ = 1.1. 
The solid line is the eigenenergy obtained from the dispersion equation (fT3|) supplemented 
by the boundary condition [Eq. ( l22l) ]. The dashed line in FIG. [2] was obtained from the 
EWS definition of the scaling energy, E sc , of Eq. (TT7I) . A significant upward shift of the 
exact eigenenergy (solid line) with respect to the scaling one, E sc , is due to the Fermi surface 
distortions of the higher multipolarities (I > 2) contributed to the dispersion equation ffTBl . 
The values of the Landau amplitude F[ = 1.1 and the parameter Q = 10.5 MeV have been 
chosen such that the best fit of the A dependence of the eigenenergies to the experimental 
data is obtained. 

Let us now consider the enhancement factor n of the isovector EWS mi. For infinite 
nuclear matter, it is given by the value of knm in Eq. ( 1T61) . The experimental determination 
of the enhancement factor k is connected to the investigation of the photoabsorption cross 
section cr a bs(^) of 7 quanta. For the velocity- independent forces, the isovector EWS m\ is 
model independent and reads (TRK sum rule) [3] 

00 

f „_ , . , 27i 2 he 2 NZ . . 

m^TRK = / d(hW) OabsM = —. (23) 

J mc A 



The photoabsorption cross section cr a b s (cf) can be expressed in terms of the strength function 
S(u) = lmx(uj)/iT as follows [21] : 

, . 47r 2 e 2 NZ n . . 

aabsM = mMio—" (24) 
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where the wave number qo{A) has to be found from Eq. (1221) in the limit of velocity- 
independent forces, i.e., at Fi = and F[ = 0. 

In a general case of velocity-dependent nuclear forces, by using Eq. (TTo^) for m 1 and Eqs. 
and ([21]), we generalize the TRK sum rule in the form 



mi(A) 



27r 2 he 2 NZ 



mc 



A 



qo(A)\ 



[1 + knm) 



(25) 



where q x (A) is derived by Eq. ([22]) at F x ^ and F[ ^ 0. 

In FIG. El the EWS enhancement factor 1 + k(A) is plotted as a function of the mass 
number A: 



mi(A) 



qo(A) 



[1 + k nm ) = 1 + k(A). (26) 

The exceedance of the 100% sum rule for rhi(A), which is observed for the IVGDR, is caused 
by the dependence of the effective nucleon-nucleon interaction on the nucleon velocity. For 
the value of the isovector amplitude F[ = 1.1, one can fit (on average) the results of our 
calculations of 1 + k{A) (solid curve in FIG. [3]) to the experimental data of the Livermore 



group 



22] . Our estimate to the enhancement factor k(A) is about 10% for light nuclei and 



increases to 20% for heavy nuclei. This result is consistent with the general conclusion of 
Ref. 22] ; that is, for A > 100, the experimental data center is at about 1.2 TRK sum- 
rule units. For the mass region A < 70, the TRK sum rule is not exhausted. As noted 
in Ref. 22], this no doubt results from the neglect of the (7,p) channel contribution for 



these nuclei. It was reported recently that inclusion of the contribution from the (7,p) cross 
section increases the EWS exhaustion from 0.87 to 1.15 for 60 Ni and from 0.64 to 0.92 for 
63 Cu 23]. The corresponding new data are shown in FIG. [3] by the symbol *. The result 



of a microscopic R 



'A calculation of the enhancement factor is shown in FIG. [3] as a dashed 



line (see Refs. 00). We point out that both our Fermi-liquid approach and the RPA 
prediction give very similar A dependence of the enhancement factor 1 + n(A). However, it 
is seen from FIG. |3]that the RPA calculation overestimates the magnitude of k{A) (see also 
Refs. (3), 8, 2I). The RPA result for k(A) can be improved by a fit of the relevant parameters 
t\ and ti in the Skyrme forces. 

The corresponding study has been recently performed within the microscopic Hartree- 
Fock plus RPA approach in Ref. j^, where a best value for n(A) of 0.22 ± 0.04 was deduced 
from a fit to the experimental data for the nucleus 208 Pb. However, reducing the enhance- 
ment factor by the variation of the Skyrme force parameters, one can expect a significant 
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change of the A behavior (slope) of the corresponding curve n(A) (dashed line in FIG. [3]). 
Unfortunately, in the nuclear literature, this kind of analysis has not yet been carried out. 
We also note that the latest microscopic RPA calculations presented in Refs. {4, 5] are 
related to nuclear matter and the nucleus 208 Pb only. They do not show an A dependence 
of the enhancement factor, but they do give a demonstration of significant variation of the 
dependence of the enhancement factor on the choice of the Skyrme force parametrization 
(see the last column in Table I of Ref. [51]). 

V. CONCLUSIONS 

In conclusion, we wish to comment that it was conceptually important for us to achieve 
a description of the A dependence of both the IVGDR energy and the enhancement factor 
simultaneously (see FIGs. 2 and 3). In our approach, we used appropriate boundary condi- 
tions that allowed us to combine both the Steinwedel- Jensen and Goldhaber- Teller models. 



A similar problem was considered earlier in Ref. [13(] but within the scaling approximation 
only. Thus, an inclusion of the effective isovector surface stiffness Q into the boundary con- 
dition [Eq. fl22|) ] leads to the A dependence of the value qRo, which becomes significantly 
smaller than the Steinwedel- Jensen's estimate qR = 2.08 p. Using the obtained value of 
qRo, we have described the IVGDR energies within the Landau kinetic theory quite well. 
Fitting the slope of the energy dependence on the mass number A to the experimental data, 
we have estimated the value of the effective isovector surface stiffness as Q ~ 11 MeV. 

The dependence of the effective nucleon-nucleon interaction on the nucleon velocity causes 
the 100% exhaustion of the TRK sum rule to be exceeded for the IVGDR. Within Landau 
kinetic theory, the EWS enhancement factor k^m m infinite nuclear matter depends on 
the interaction amplitudes Fi and F[. In finite nuclei, the A dependence of the EWS 
enhancement factor k(A) occurs because of the boundary condition [Eq. ( |22i) ]. The value of 
k(A) increases with A. A fit of the enhancement factor to the proper experimental data leads 
to a value for the isovector Landau amplitude of F[ ~ 1.1. The obtained value of F[ exceeds 
;he estimate F[ = 0.5 — 0.7 derived earlier from Skyrme forces for infinite nuclear matter 
9, [ll|, [if]]. This exceedance appears since our derivation of F[ is related to the interior of 



the finite nucleus. For finite nuclei, the Landau amplitudes and F' k are r-dependent ones 
with a bump within the nuclear surface fill . \i\ . This effectively increases the bulk values 
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of Ff. and F' k in the limit of the sharp nuclear surface assumed in this paper. 

We show that the value of the enhancement factor is about 10% for light nuclei and 
reaches approximately 20% for heavy nuclei. For nuclei with A > 100, our results are close 



to the experimental data from Livermore discussed in Ref. [22j and they are in agreement 
with the best value for the enhancement factor obtained for the nucleus 208 Pb in Ref. 
We note also that the A behavior of the enhancement factor obtained in our Fermi-liquid 
approach is similar to the one derived from the microscopic RPA calculations (see dashed 
line in FIG. 3), but, as was previously mentioned, the RPA results show a strong variation 
of the enhancement factor with the choice of the Skyrme force parametrization [5(. 
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FIG. 2: Dependence of the IVGDR energy on the mass number A. The dashed line is the calculation 
that includes Fermi surface deformations up to quadrupole order [scaling approximation: see E sc 
in Eq. (|17|) ]; the solid line was obtained from the dispersion equation and the secular equation 



22]) . The dots are the experimental data taken from Ref. 



221. 
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FIG. 3: Dependence of the enhancement factor 1 + k(A) for IVGDR on the mass number A. The 
calculation result was obtained for a Landau amplitude F\ = —0.3, F[ = 1.1 (solid curve) . The 

nn 

dashed line is the microscopic RPA calculations with Skyrme forces SkM* from Refs. [1], HQ]. The 
solid points are the experimental data of the Livermore group from Ref. [22] . Two points (noted 
by the symbol *) were obtained by the inclusion of the contribution from the (j,p) cross section. 
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